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Abstract 

We perform the analysis of a hyperbolic model which is the analog of the Fisher-KPP 
equation. This model accounts for particles that move at maximal speed (e > 0), and 
proliferate according to a reaction term of monostable type. We study the existence and 
stability of traveling fronts. We exhibit a transition depending on the parameter e: for small 
e the behaviour is essentially the same as for the diffusive Fisher-KPP equation. However, 
for large e the traveling front with minimal speed is discontinuous and travels at the maximal 
speed e~^. The traveling fronts with minimal speed are linearly stable in weighted spaces. 
We also prove local nonlinear stability of the traveling front with minimal speed when e is 
smaller than the transition parameter. 

Key-words: traveling waves; Fisher-KPP equation; telegraph equation; nonlinear stability. 



1 Introduction 

We consider the problem of traveling fronts driven by growth (e.g. cell division) together with cell 
dispersal, where the motion process is given by a hyperbolic equation. This is motivated by the 
occurence of traveling pulses in populations of bacteria swimming inside a narrow channel [TJ . 
It has been demonstrated that kinetic models are well adapted to this problem j37j . We will focus 
on the following model introduced by Dunbar and Othmer [T^] (see also Hadeler Plj) and Fedotov 

[Biiaiis] 

e^dttPs{t,x) + {l~e^F'{p,it,x)))dtp,{t,x)-d,.,p,{t,x)^F{p,{t,x)), t>Q, xGR. (1.1) 

The cell density is denoted by Pe{t, x). The parameter e > is a scaling factor. It accounts for the 
ratio between the mean free path of cells and the space scale. The growth function F is subject 
to the following assumptions (the so-called monostable nonlinearity) 

F is concave , F{0) = F(l) = , F{p) > if p G (0, 1) . 

For the sake of clarity we will sometimes take as an example the logistic growth function F{p) — 

pa-p)-^ 

Equation (1.1) is equivalent to the hyperbolic system 

dtPe + e-^d^ (je) = F{p^) 
edtje + d^pe = -e^^je . 



We will successively use the formulation (1.1 1 or the equivalent formulation (1.2) in the sequel. 
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Since the pioneering work by Fisher |19j and Kolmogorov-Petrovskii-Piskunov [27], dispersion 
of biological species has been usually modelled by mean of reaction-diffusion equations. The 
main drawback of these models is that they allow infinite speed of propagation. This is clearly 
irrelevant for biological species. Several modifications have been proposed to circumvent this issue. 
It has been proposed to replace the linear diffusion by a nonlinear diffusion of porous-medium type 
pS8 ,.29, 33 . This is known to yield propagation of the support at finite speed [30, 31 . The density- 
dependent diffusion coefRcient stems for a pressure effect among individuals which influences the 
speed of diffusion. Pressure is very low when the population is sparse, whereas it has a strong 
effect when the population is highly densified. Recently, this approach has been developped for 
the invasion of glioma cells in the brain [S]. Alternatively, some authors have proposed to impose 
a limiting flux for which the nonlinearity involves the gradient of the concentration IH [7l|4] . 

The diffusion approximation is generally acceptable in ecological problems where space and 
time scales are large enough. However, kinetic equations have emerged recently to model self- 
organization in bacterial population at smaller scales [3J [351 [T31 [IHl [32 [Ml ISZ] • These models are 
based on velocity-jump processes. It is now standard to perform a drift-diffusion limit to recover 
classical reaction-diffusion equations [25l [9j [131 [22] . However it is claimed in [37] that the diffusion 
approximation is not suitable, and the full kinetic equation has to be handled with. 

Hyperbolic models coupled with growth have already been studied in [T^ [211 [211 [ID] ■ In [21] it 
is required that the nonlinear function in front of the time first derivative dtPe is positive (namely 
here, 1 — e^F'{u) > 0). Indeed, this enables to perform a suitable change of variables in order to 
reduce to the classical Fisher-KPP problem. In our context this is equivalent to e^F'(O) < 1 since 
F is concave. In [21j this nonlinear contribution is replaced by 1, the authors study the following 
equation (damped hyperbolic Fisher-KPP equation): 

£^dttPe{t,x) + dtPe{t,x) - dxxPE{t,x) = F{pe{t,x)) . 



We also refer to [ID] where the authors analyse some kinetic equation more general than ( 1.1 1 
in the regime of small parameter e, i.e. a perturbation of the diffusive regime. 

It is worth recalling some basic results related to reaction-diffusion equations. First, as e — 



the density solution to (1.1) formally converges to a solution of the Fisher-KPP equation: 

dtPQ{t,x) - dxxPo{t,x) = F{po{t,x)) . 

The long time behaviour of such equation is well understood since the pioneering works by 
Kolmogorov-Petrovsky-Piskunov [27^ and Aronson- Weinberger [5]. For nonincreasing initial data 
with sufficient decay at infinity the solution behaves asymptotically as a traveling front moving 
at the speed s — 2y/F'{0). Moreover the traveling front solution with minimal speed is stable in 
some weighted space [20]. 

In this work we prove that analogous results hold true in the parabolic regime e'^F'{0) < 1. 



Namely there exists a continuum of speeds [s*(£),e ^) for which (1.1 1 admits smooth traveling 
fronts. The minimal speed is given by |14j 



Obviously we have s*(e) < min(2-\/i^'(0), e ^). There also exists supersonic traveling fronts, with 
speed s > This appears surprising at first glance since the speed of propagation for the 



hyperbolic equation ( 1.1 1 is e ^ (see formulation (1.2) and Section^. These fronts are essentially 



pulled by growth, since they travel faster than the maximum speed of propagation. The results 
are summarized in the following Theorem. 

Theorem 1 (Parabolic regime). Assume thate'^F'{0) < 1. Then for all s <E [s* (e), e~^) , there 
exists a smooth traveling front solution of ( 1.1 ) with speed s. For all s G (e^^, oo) there also exists 



a smooth traveling front. For s — e there exists a weak traveling front. 
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We also prove that the minimal speed traveling front is nonlinearly stable (see Section [sj 
Theorem 13 ). 

There is a transition occuring when e^i^'(O) = 1. In the hyperbolic regime £^F'{0) > 1 the 
minimal speed speed becomes: 

s*{e)^e-\ if£2F'(0)>l. (1.4) 

On the other hand, the front traveling with minimal speed s* (e) is discontinuous as soon as 
e^F'(O) > 1. In the critical case £^F'(0) = 1 there exists a continuous but not smooth traveling 
front with minimal speed s* — \J F'{0). 

Theorem 2 (Hyperbolic regime). Assume that e^F'(O) > 1. Then there exists a weak traveling 
front solution of ( |1.1[ ) with speed s*{e) — . The wave is discontinuous if e^F'{0) > 1. For all 
s € (e^^jOo) there exists a smooth traveling front. 



We conclude this introduction by giving the precise definition of traveling fronts (smooth and 
weak) that will be used throughout the paper. 

Definition 3. We say that a function p{t,x) is a smooth traveling front solution with speed s 
of equation (1.1 1 if it can be written p{t,x) = v{x — st), where v G C^(M), > 0, i^(— oo) = 1, 



;/(+cx)) = and v satisfies 

(e^s^ - l)z/"(z) - (l - e2F'(iy(z))) siy'(z) = F(i/(z)) , zeR. (1.5) 
We say that p is a weak traveling front with speed s if it can be written p{t, x) ~ v[x — st), where 



V G i°°(IR), v>Q, v{—oo) = 1, i^(+oo) = and v satisfies (1.5) in the sense of distributions: 

-F{v)Lp\dx = 0. 



V(p e V{M) I ((e^s^ - l)vip" +{u^ se^F{u)) sip' - 



In the following Section[2]we show some numerical simulations in order to illustrate our results. 
Section [S] is devoted to the proof of existence of the travling fronts in the various regimes (resp. 
parabolic, hyperbolic, and supersonic). Finally, in Section |4] and Section [5] we prove the stability 
of the traveling fronts having minimal speed s*{e). We begin with linear stability (Section |4]) since 
it is technically better tractable, and it let us discuss the case of the hyperbolic regime. We prove 
the full nonlinear stability in the range e S (0, 1/ y^F'{0)) (parabolic regime) in Section [i] 

2 Numerical simulations 



In this Section we perform numerical simulations of (1.1). We choose a logistic reaction term: 
F{p) = p{l — p). We first symmetrize the hyperbolic system (1.2) by introducing /+ — \{p + j) 
and = \{P ^ j)- This results in the following system: 

dj+{t,x)+e-^dj+it,x) = ^ {f-{t,x) - f+it,x)) + lF{p{t,x)) 

: (2.1) 

dtf-it,x) - e-^d,f'{t,x) = if+it,x) - f-{t,x)) + -F{pit,x)) . 

In other words, the population is split into two subpopulations: p = f^ + where the density 
/+ denotes particles moving to the right with velocity e~^, whereas /~ denotes particles moving 
to the left with the opposite velocity. 

We discretize the transport part using a finite volume scheme. Since we want to catch discon- 
tinuous fronts in the hyperbolic regime e^F' {0) > 1, we aim to avoid numerical diffusion. Therefore 
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we use a nonlinear fiux-limiter scheme [231 HI] • The reaction part is discretized following the Euler 
explicit method. 

Jn + \,i ^ Jn,i ^ \Jn,i^Pi 2 ■'n,i-l A-1 2 / 2 ^ 2 ' 

The non-linear reconstruction of the slope is given by 

-, I fn,i ^ fn,i-i -fn.i+i^fti] ■ ^ jO if sign (p) 7^ sign (9) 

Pt = mmmod , ^ . , mmmod (p, (?) = < . /, , , . / ^ ... , . ■ , . 

\ /\x iXx J I mm(|p|, |q|)sign (p) if sign (p) = sign (9) 

We compute the solution on the interval (a, b) with the following boundary conditions: f^(a) = 1/2 
and f~{b) = 0. The discretization of the second equation for /~ (2.1) is similar. The CFL 
condition reads At < eAx. It degenerates when e \, 0, but we are mainly interested in the 
hyperbolic regime when e is large enough. Other strategies should be used in the diffusive regime 
£ ^ 1, e.g. asymptotic-preserving schemes (see [IZIIH] and references therein). 



3 Traveling wave solutions 
3.1 Characteristic equation 



We begin with a careful study of the linearization of ( 1.5 ) around ^ 0. We expect an exponential 



decay e as z ^ +00. The characteristic equation reads as follows, 

{eh^ - 1)A2 + (1 - e^F'{0))sX - F'(0) = . (3.1) 

The discriminant is A = (e'^F'(O) + l)^ s'^ — 4F'(0). Hence we expect an oscillatory behaviour in 
the case A < 0, i.e. s < s*{e). We assume henceforth s > s*{e). In the case s < {subsonic 
fronts) we have to distinguish between the parabolic regime e'^F'{0) < 1 and the hyperbolic regime 



e F'(0) > 1. In the former regime equation (3.11 possesses two positive roots, accounting for a 



damped behaviour. In the latter regime equation (3.1) possesses two negative roots. In the case 
s > {supersonic fronts) we get two roots having opposite signs. 

Next we investigate the linear behaviour close to v 1. We expect an exponential relaxation 
1 — ^ as z — > —00. The characteristic equation reads as follows, 

{e^s^ - 1)A'2 - (1 - s^F'{l))sX' - F'{1) = . (3.2) 



We have A' = [e^F'{l) + IJ^s^ - 4F'(1) > 0. In the case s < e^^ equation ( |3.2| possesses two 
roots having opposite signs. In the case s > e^^ it has two positive roots. 

We summarize our expectations about the possible existence of positive traveling fronts in 
Tabled 

3.2 Obstruction for s < s*{e) 

In this section we prove that no traveling front solution exists if the speed is below the minimal 
speed. 

Proposition 4. There exists no traveling front with speed s for s < s*{e), where s*{e) is given 



by (1.3)-(1.4) 



Proof. We argue by contradiction. The obstruction comes from the exponential decay at -t-00. 
Assume that there exists such a traveling front ^{z). Classical regularity estimates show that v 
is smooth. It is necessarily decreasing as soon at it is below 1. Otherwise, it would reach a local 
minimum at some point zq G K, for which v'{zo) = and F{u{zo)) < 0. As iy{zo) < 1, this would 
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Figure 1: Numerical simulations of the equation (1.1 ) for F{p) = p{l — p) and for different values 
of £ = 0.5, 1, 2. Numerical method is described in Section[2] The initial data is the step function 
f^{x < 0) = 1, f^{x > 0) = 0, and /~ = 0. For each value of e we plot the density function 
p = + ,f^ in the (x, t) space, and the density p{tQ, •) at some chosen time to. We clearly observe 
in every cases a front traveling asymptotically at speed s*(e) as expected. We also observe the 
transition between a smooth front and a discontinuous one. The transition occurs at e = 1. In the 
case e = 1 we have superposed the expected profile (l — e^^^)_|^ in black, continous line, 

for the sake of comparison. 
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parabolic 



s < e 



if s < s*(e), NO 



if s > s*(e), YES 
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s > e~ 



YES 




hyperbolic 




YES 




Table 1: Phase plane dynamics depending on the regime (parabolic vs. hyperbolic) and the value 
of the speed with respect to s*{e) and e~^. Note that in the supersonic case we have reversed 
the sign of the arrows to clarify the picture. In every picture the red line represents the traveling 
front trajectory, and the green lines are the axes {u — 0} and {v = 0}. We do not consider the 
case s = since the dynamics are singular in this case and should be considered separately (see 



Section 3.4 1 



give v^zq) = 0. Uniqueness would imply 1^ = 0. Next, we define the exponential rate of decay at 
+oo: 

A := liminf — -V > 0. 
Consider a sequence Zn — >■ +oo such that — z^'(z„)/i^(z„) A and define the renormalized shift: 



iy„{z) 



v{Zn) 



This function is locally bounded by classical Harnack estimates. It satisfies 

[e^s^ - \)vl{z) + {e'F\v{z + z„)) - l) sv'^z) = -^F{v{zn)Vn{z)) , z G M . 

V\Zn) 

Up to extraction, we can assume that {h'n)n converges to some function Vrx, in ^loci^)- T^^^ limiting 
function satisfies 



{e's' - iy:,{z) + {e'F'iO) - l) sj^U^) = F'{0),y^{z) , z € 



(3.3) 



Define V = v'^/voo- The definition of A yields minR V = F(0) = -A. Thus V'{Q) = 0. Hence we 
deduce from (3.3| that v^oiz) = Voo{^^)e^^^' ■ Plugging this into (3.3), we obtain that A satisfies 
the following second order equation. 



We know from Section [3J] that there is no real root in the case s < s*{e). 



□ 
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3.3 Existence of smooth traveling fronts in the parabolic regime s G 

In the author prove the existence of traveling front, by reducing the problem to the classical 
Fisher-KPP problem. It is required that the nonlinear function l — e'^F'{u) remains positive, which 
reads exactly e^F'(O) < 1 in our context. We present below a direct proof based on the method 
of sub- and supersolutions. 

3.3.1 The linearized problem 



Proposition 5. Let \g be the smallest (pos itive ) root of the characteristic polynomial (3.1). Then 
^{z) = min{l, e""**"^} is a supersolution of (1.5). 

Proof. Let r(z) — e'~^'' . Then as r is decreasing and F is concave, it is easy to see that r is a 
supersolution of (1.5). On the other hand, the constant function 1 is clearly a solution of (1.5). It 
is clear that the minimum of two supersolutions is a supersolution. □ 

3.3.2 Resolution of the problem on a bounded interval 
Proposition 6. For all a > and r G M, there exists a solution Va.r of 

'^a^rl-o) = I^(-a + r), (3.4) 
Va.ria) ^ V{a + t). 

Moreover, this function is nonincreasing over {—a, a) and it is unique in the class of nonincreasing 
functions. 

In order to prove this result, we consider the following sequence of problems: 

• vq{z) ~ V{x + r) 

• Vn+i is solution to 

■(£2s2 - l)i.;'^^ + {e'^F'iun) - l)sj/;+i + Mv^+i = F{u„) + Mvn in (-a, a), 
M„+i(-a) = I7(-a + t), (3.5) 
Vn+i{a) = v{a + t), 

where M > is large enough so that s i~-> F{s) + Ms is increasing. 

Lemma 7. The sequence {yn)n is well-defined. The functions z i— !■ t'„(z) are nonincreasing and 
for all z e (—a, a), the sequence (j^n(z))n is nonincreasing. 

Proof. We prove this Lemma by induction. Clearly, i^q is nonincreasing. Let wq = I'l — vq. As vq 
is a supersolution of equation (1.5), one has 

(£2^2 _ {£^F'{va) - l)sw'a + Mwq < in (-a, a), 

wo{-a) = wo{a) = 0. 

As M > 0, the weak maximum principle gives wq < 0, that is, j^i < vq. 

Assume that Lemma [7] is true up to rank n. Let Wn = I'n+i — t^n- As F is concave and 
Vn-i ^ t^ni "we know that i^'(j/„_i) < F'{vn). As Vn is nonincreasing, we thus get 

{e^s^ - IX + {e^F'{v,,) - l)s< + Mwn < in (-a, a), 
Wn{-a) = Wn{a) = 0. 

Hence, Wn < and thus I'n+i < i^n- 
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Next, as F{v{a + t)) > 0, the constant function v{a + t) is a subsolution of (3.5 1. Thus the 
same arguments as before give 

V{a + r) < Vn+i{z) < vq^z) = v{z + t) < v{—a + t) in [—a, a] 

and the strong maximum principle yields that these inequalities are strict in (—a, a). 

Fix A > and define Sa := (a — -^j a) and := i^„+i(A — 2a+z). The previous inequalities 

yield that Vn+i < i^n+i fo^' small A > 0. Let now increase A and set 

A* = sup{A e (0, 2a), i^n+i < i^n+i in for M ^ (0, A)} > 0. 

Assume that A* < 2a. Then classical continuity arguments yield that there exists z_ £ Sa» such 
that Vn+i{z) = v^j^i{z). Define v = ly^+i ~ i^n+i- Then this function is nonnegative, vanishes in 
z and, as F is concave and u^+i is nonincreasing with respect to z, it satisfies 

(e2s2 - l)v" + {e'^F'{un+i) - l)sv' + Mv < F{i^^') - F{vr,) + M(z/f - z/„) < in Sa-, 
v[a - A*) = Vn+ii-a) - Un+iia - A*) > 0, 
v[a) = ^'„+i(A* - a) - Vn+i{a) > 0. 

Hence, the strong maximum principle gives that v = v{z) = 0, which contradicts the boundary 
conditions. Thus, A* = 2a and thus < i^n+i in for all fi G (0, 2a), which means that Vn+i 
is nonincreasing. □ 

Proof of Proposition^ As the sequence (vn)n is decreasing and bounded from below, it admits a 
limit Va.T as n -> +oo. It easily follows from the classical regularity estimates that Va^r satisfies 
the properties of Proposition |6] 

If vi and V2 are two nondecreasing solutions of (3.4), then the same arguments as before give 
that i/f < vi in for all /i e (0, 2a). Hence, vi < 1/2 and a symmetry argument gives 1/1 = 1/2- Q 

Lemma 8. For all a> 0, there exists e M such that i'a,T„(0) = \- 

Proof. Define I{t) := i^a,r(0). Then it easily follows from the classical regularity estimates and 
from the uniqueness of Va^r that / is a continuous function. Moreover, as i'a,T is nonincreasing, 
one has 

V{a + t) < I{t) < v{-a + t). 

As + r) — > as T — > +00 and + r) — > 1 as r — >■ —00 locally uniformly on M, one has 
/(— 00) = 1 and /(+00) = 0. The conclusion follows. □ 

3.3.3 Existence of traveling fronts with speeds s E [s*{e),e^^) 

We conclude by giving the proof of Theorem [T] as a combination of the above results. 

Proof of Theorem^ The classical elliptic regularity estimates yield that there exists a sequence 
a„ — !■ +00 such that (i'a„,r„ )n converges in C;2g^(M) as n — > +00 to a function v. Then 1/ satisfies 
(1.5), it is nonincreasing, < < 1 and 1/(0) = 1/2. It is easy to see that 00) = 1 and 
i/{+oo)=0. □ 

We will need the following classical inequality satisfied by the traveling profile. 

Lemma 9. The traveling profile v satisfies: Vz v' {z) + Xv{z) > 0, where A is the smallest positive 
root of (3.1 ). 

Proof. We introduce (p{z) = —77^. It is nonnegative, and it satisfies the following first-order 
ODE with a source term 



{s's' - 1) (-^'(z) + ^{zf) + (1 - s'F'ii/{z))) s^iz) = 



iy{z) 
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Since F is concave, ip satisfies the differential inequality 

(1 - e^s^) (^'(z) < (1 - £2s2) <^(z)2 _ (1 _ £2^'(0)) + F'(0) . 



The right-hand-side is the characteristic polynomial of the linearized equation (3.1 ). Moreover the 
function verifies Ivcuz^-ao ^(z) = 0. Hence a simple ODE argument shows that Vz ip{z) < X. □ 

3.4 Existence of weak traveling fronts in the hyperbolic regime s = e^^ 

In this Section we investigate the existence of traveling fronts with critical speed s = s^^ in the 
hyperbolic regime £^F'(0) — 1. 

Proof of Theorem^ We distinguish between the two cases £^F'(0) > 1 and e^i^'(O) = 1. 
1- First case: e'^F'{0) > 1 

Using the concavity and continuity properties of F, we define G (0, 1) as follows: 

6=: is the unique positive root of e'^F{p) = p. 
By concavity of F we get that e'^F'(p) < 1 for all p> 9e- Set v the maximal solution of 

eF{v{z)) 




e'^F'{v{z)) - 1' 



(3.6) 



2 

Let / = (— oo, zi) be the (maximal) interval of definition of v, with G /. Let 

zq = sup{z e /, v{z) > 9^}. 

Hence we have i^iz) — ^ 0^ as z -->■ Zq . Moreover, as v is clearly decreasing on (—00,0), one easily 
gets i'{—oo) = 1. 

We set v{zq) — 9^ and we extend by over (zq, oo). Then it is easy to see that v is a weak 
solution, in the sense of distributions, of 



{e^F{iy)-iyy = sF{v) 



on 



since £^^(0) = and e^F{9,) = 9^. 

Up to space shifting z — zq, we assume that the discontinuity arises at z = 0. 

Exemple: the case F{p) = p{l — p) and £ > 1. The traveling profile solves 

^ ev{z){l - v{z)) 

^ ' £2 - 1 - 2£2z.(z) ' 

or equivalently 

v{zy^-^ (1 - iyiz)f+^ = ke" . 

The constant k is determined by the condition iy{0) 6'^ = 1 — £^^. Finally the traveling profile 
v{z) satisfies the following implicit relation: 

,y{zf-^{l-l^{z))'+'" = (l-£-2)^'"l(£-2)-' + le- = (£2-l)^'-^"+2-'l°g^\ (3.7) 
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2- Second case: e^F'{0) = 1. 

It is easy to see that e (l/i^'fO), +00) 1— > 6^ is continuous and nonincreasing. Hence, it admits 
a limit 6 when e — > 1/^F'(0)(N). We consider the family of weak traveling front solutions of 
). Then, the family {i'e)e converges as e — >■ 1/-^/ -^'(0) to a measurable nonincreasing function 
lich satisfies 

1 \' 1 

F{iy) -1^] = F{v) on M 



(3.6 



,^'(0) ' ' J /P(0) 
in the sense of distributions. This is the desired profile 



□ 



Example: the case F[p) = p(l — p) and e = 1. The traveling profile reads (3.7): 

K.) = (l-e//^)^. 

3.5 Existence of supersonic traveling fronts s > 

In this Section we investigate the existence of supersonic traveling fronts with speeds above the 
maximal speed of propagation s > e~^. These fronts are essentially driven by growth. The 
existence of such " unrealistic" fronts is motivated by the extreme case e — > +00 for which we have 
formally dtp — F{p) \1.2\ . There exist traveling fronts of arbitrary speed which are solutions to 
-sv' = F{v). 

Proposition 10. Given any speed s > e^^ there exists a smooth traveling front v{x — st) with 
this speed. 

Proof. We sketch the proof. We give below the key arguments derived from phase plane analysis. 
The same procedure as developped in Section |3.3| based on sub- and supersolutions could be 
reproduced based on the following ingredients. 



We learn from simple phase plane considerations associated to (1.51 that the situation is re- 
versed in comparison to the classical Fisher-KPP case (or £^F'(0) < 1 and s G [s* (e) , e~^)) . 
Namely the point (0,0) is a saddle point (instead of a stable node) whereas (1,0) is an unstable 
node (instead of saddle point). This motivates "time reversal": V{z) — i>{~z). Equation ( |1.5[ ) 
becomes 



l)V"{z)-{e^F'{Viz))~l)sV\z)=F{Viz)) , ze 



[e s 



We make the classical phase-plane transformation V' = P [27l [18] . We end up with the implicit 
ODE with Dirichlet boundary conditions for P: 

ie's' - \)P\v) - {e'F' („) - l) . = ^ , P(0) = P(l) = . 

The unstable direction is given by P{v) — Xv where A is the positive root of 

{eh^-l)X-{e^F'iO)-l)s^^. (3.8) 

A 

Since F is concave we deduce that P{v) = Au is a supersolution as in Proposition [sj In fact, 
denoting Q{v) = P{v) — Xv we have 

(£2^2 _ ^)g,(^) ^ (^2^2 _ l)(p'(^,) _ At-) < F'{0)V ( - \v) < i?'(0) -"^^"^ 



P{v) X J - ' ' P{v)X 

Hence the trajectory leaving the saddle point (0, 0) in the phase plane (V, V') remains below the 
line V' < XV. On the other hand it is straightforward to check that kF{v) is a supersolution 



^If F{s) < F'{0)s when s > is small, then one can prove that 8 = 0. But this is not true anymore if F is linear 
in the neighbourhood of 0. 
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Figure 2: Supersonic traveling front in the phase plane (V, V) for the nonlinearity F{p) = p{l — p), 
and parameters e — \/2, s = 1 > Be aware of the time reversal iy{z) — V{—z), which is the 
reason why V' > 0. The red line represents the traveling profile, and the green line represents the 
supersolution kF{v). 



where k = e^s/(e^s^ — 1). We denote R{v) = P{v) — kF{v). We have ks > 1 and 

(£^52 _ = (^2^2 _ _ kF'{v)) ^ ehF'iv) - s + ^ - (£^52 _ i)kF'{v) 

P{v) 

1 R(v) R(v) 
= -s + - - , „\ \ < ^ ' 



k kP{v) kP{v) 

We also show that initially (as w — > 0) we have kF'{0) > A. This proves that R{v) < for all 
V G (0, 1). Indeed, we plug fcF'(O) in place of A into (3.8 1 and we get 

ie^s^ - l)fc^^'(O) - (£^F' iO)-l)s = s>l^^^. 

As a conclusion the trajectory leaving the saddle node at (0, 0) is trapped in the set {0 < ?; < 
1 1 < p < kF{v)} (see Fig. [2]). By the Poincare-Bendixon theorem it necessarily converges to the 
stable node at (1, 0). This heteroclinic trajectory is the traveling front in the supersonic case. □ 



4 Linear stability of traveling front solutions 

In this Section we investigate the linear stability of the traveling front having minimal speed 
s ~ s*{e) in both the parabolic and the hyperbolic regime. We seek stability in some weighted 
L2 space. The important matter here is to identify the weight e"^. The same weight shall be used 
crucially for the nonlinear stability analysis (Section [5]). 
We recall that the minimal speed is given by 

1^4^ if£^nO)<l 
s (e) = \ 1 + £2F'(0) ^ ' 

if£2^^'(0)>l 

The profile of the wave has the following properties in the case £2i^'(0) < 1: 
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where the decay exponent A is 

_ s(l - e^f ^(0)) _ l + e^F'{0) 
^ 2(1 - £2.s2) - 1 - £2^'(0) ■ 

We will use in this Section the formulation ( |1.2| of our system. The linearized system around 
the stationary profile v in the moving frame z — x — st reads 

{dt-sd,)u + dj-) ^F'{v)u 

„ (4.1) 
e{dt - sdz)v + dxu = . 

e 

Proposition 11. Let e > 0. In the case > 1/F'{0) assume that the initial perturbation has the 
same support as the profile. There exists a function (p^iz) such that the minimal speed traveling 
front is linearly stable in the weighted L'^{e'^'^^^^^dz) space. More precisely the following Lyapunov 
identity holds true for solutions of the linear system ( |4.1[ ) . 

d_ /I 
Jt \2 

Proof. We denote (j) = (j)^ for the sake of clarity. We multiply the first equation by ue'^'^, and the 
second equation by ve'^'^, where (j) is to be determined. We get 



{\u\^ + \v\^)e^*^^'Uz) <0. 



d /I 
dt 

d /I 
dt \ 2 



(I I |?^pe2^(^)dz) + J / \u\^d,e^'^'^'Uz+ [ dJ-)ue^'f'^'Uz= [ F' iiy)\u\^ e^'f'^'Uz , 
V 2 Jr / 2 Jr Jr ^ £ ^ Jr 

vl^e^^^'Uz^ + J \v\'d,e^*<^'Uz + J ve^'>''^'Uz = j 



Summing the two estimates we obtain 

+ / {sd,^{z)-F\v))\u\''e^'^'^'Uz+ I {sd,^{z) + ^ \v\'' e'^^'Uz-- [ {d,(l){z)) uve^^^'-^z = . 
Jr Jr\ £ / ^ Jr 

We require that the dissipation is nonnegative. This is guaranteed if the following discriminant is 
nonpositive: 

A(z) = 4 {d^m f ~ 4 {sd,cj){z) - F'i,,)) (sd^z) + ^ 



((1 - e's') {dMz)f - ,s (1 - e^F'iv)) d,dp{z) + F'{v)) . (4.2) 



The rest of the proof is devoted to finding such a weight (j){z) satisfying this sign condition. We 
distinguish between the parabolic and the hyperbolic regime. 

1- The parabolic regime. 

In the case £^^"(0) < 1 we have e^s"^ < 1. Hence the most favourable choice for the discriminant 
(4.2) to be negative is: 

s{l-e^F'{u)) ^ 1 - e^F'jv) ^ ^ 
2(1 -£2.2) = \-s2F^(0) ^0- 
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Notice that dz(p — > A as z — > +00 [21]. We check that the discriminant is indeed nonpositive: 



A(z) = -4(1- e^s^) {d,^{z)f + 4F'{iy) 



e 



(-4F'(0) (1 + e^F'{v)f + ^F'{v) (l + £^^^'(0))') 



(l-£2^'(0))' 



— {F'{Q) - F'{v)) (1 - e^F'{Q)F'{u)) . 



{l-e^F'{0)y 

We have A(z) < since Vz F'{v{z)) < F'{0) and e'^F'{0) < 1. We have henceforth 



< J (s9,0(z) - F'(zy) - ^(z)) lupe^-^W dz + J (^sd,<j>{z) + ^ - A{z)^ \v\^e^^^'Uz, 



where A{z) is solution of 

4 (sdAiz) - F'{v) - A{z)) {sd,cj>{z) + A{z) ] = ^ {d,4>{z)f 
A straightforward computation gives 



2A{z) = ( 2s9,0(z) - + \] - { { F' [v) + 1 ) + ^ (a,0(z)) 



1/2 



\2 



((1 - £2^2)2 (1 + S^F'{y)f + £2,s2 (1 _ e^F'{l,)f) 



£2 (1 — £2s2) £2 — e2g2 

/ / A 9 \ 

1-£2F'(z.) / / /I -£2i^'(0)V 4£2F'(0) 



n2 



£2(l-e2s2) 1 U 1 +£2^^/(0) y Vl-£2i?^(j^)y (1+ £2^/(0)) 

1-£2f'H (^_(^, fl-e^F'{0)\^ fl+e^F'{u)Y (l - £^^^'(0))' 



£2(l-£2s2) I I Vl +«^^-P''(0)y Vl"£^^'(i^)y (l+£2i^'(0)) 

We clearly have A{z) > since 

l + e^F'{iy) ^ l + e^F'{0) 

1 - e'^F'{v) - 1 - £2F'(0) ■ 

Finally we obtain in the case £2F'(0) < 1, 



d (\ 
(ft I 2 



{\ [ (IwP + e'^^'^ dz) + [ A{z) {\u\^ + e2*(^) dz < . (4.3) 
\2 Jr J Jr 



2- The hyperbolic regime. 

We assume for simplicity that the support of the stationary profile is Suppz^ = (— oo,0]. 
In the hyperbolic regime we have s — £~^, so the discriminant equation (4.2 1 reduces to 

A(z) = I {-s (1 - e^F'i,,)) 9,0(z) + F\,,)) . 

We naturally choose 

= l-£2f^(.) • 
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Within this choice for cj) we get, 



d 
dt 



(If {\u\^ + \v\^) e^^^'-Uz] + f A{z){e'F'{i^{z))u-vYe^^^''>dz^O, 

\2jz<o / Jz<0 



where the additional weight in the dissipation writes: 

Aiz) = ^ 



e2{l-e'^F'{iy{z))) 



In the case e^F'(O) > 1 we have 1 — e'^F'{v{z)) > on Suppj^ (see Section 3.4 1. Notice that 
the monotonicity of may change on Suppj/ since F' {v{z)) may change sign. We observe that 
A{z) is uniformly bounded from below on Supp v. 

In the transition case e^i^'(O) = 1, we have dz<i>{z) +oo as z 0^. We observe that 
A{z) — 7> +CX3 as z — 0~ also. 

□ 



Exemple: the case F(p) = pil — p), and e = 1. We can easily compute from Section 3.4 

0(z) = -^-log(l-e^/2) . 
Remark 12 (Lack of coercivity). 1- The parabohc regime. We directly observe that A{z) — > 



as z ^ +00 in the Lyapunov identity (4.3). This corresponds to the lack of coercivity of the linear 
operator. It has been clearly identified for the classical Fisher-KPP equation \2(A This lack of 
coercivity is a source of complication for the next question, i.e. nonlinear stability (see Section^. 
2- The hyperboUc regime. The situation is more degenerated here: the dissipation provides 
information about the relaxation of v towards e'^F'{v)u only. 

5 Nonlinear stability of traveling front solutions in the parabolic 
regime £^F'(0) < 1 

In this Section we investigate the stability of the traveling profile having minimal speed s — s*(e) in 
the parabolic regime. We seek stability in the energy class. Energy methods have been successfully 
applied to reaction-diffusion euations [Ml HH |351 [52] . We follow the strategy developped in [2T] 
for a simpler equation, namely the damped hyperbolic Fisher-KPP equation. 

Before stating the theorem we give some useful notations. The perturbation is denoted by 
u(t, z) = pit, z) — v{z) where z = x — st is the space variable in the moving frame. We also need 



some weighted perurbation w = e^'w, where is an explicit weight to be precised later (5.10). 

Theorem 13. For alle e ^0, l/-\/_F'(0)^ there exists a constant c{e) such that the following claim 
holds true: let u'^ be any compactly supported initial perturbation which satisfies 

\\dznY2 + I dz + Wd.wY + \\wY2 < c{e) , 

J Z<Zq 

then it remains uniformly small for all time < > 0; 

sup {\\dMt)\\l + / I^WP dz + WdMml + W^ml] < c{e) , 

t>0 V Jz<zo J 

and the perturbation is globally decaying in the following sense: 

(\\dzu\\l+ f \u\^dz+\\dM\l+ f e-'^^''>\w\^dz) eL^{0,+oo). 

\ J Z<Zo J Z>Zo J 
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Remark 14. 1. The additional weight e in the last contribution (weighted L? space) is 

specific to the lack of coercivity in the energy estimates. 

2. The constant c{e) that we obtain degenerates as e — > 1, due to the transition from a parabolic 
to an hyperbolic regime. 

3. We restrict ourselves to compactly supported initial perturbations to justify all integration 
by parts. Indeed the solution u(t, z) remains compactly supported for allt>0 because of the 
finite speed of propagation. The result is the same if we assume that u*^ decays sufficiently 
fast at infinity. 

Proof. We proceed in several steps. 

1- Derivation of the energy estimates. The equation satisfied by the perturbation u writes 

{dttu - 2sdtzU + s^dzzu) + (l - e^F'{v + u)) {dtu - sdzu) - dzzU 

+ £^{F'{v + m) - F'{v))sdzv = F{v + u) - F{v) . 

We write the nonhnearities as follows: 

F'{v + u) = F'{iy) + Ki{z;u)u, 

F'{v + u)- F'{v) = F"{v)u + K2{z] u)u^ , 

F{v + u)- F{v) = F'{v)u + K^iz; u)u'^ . 

where the functions Ki are uniformly bounded in i°°(M). More precisely we have 

Ki{z;u)^ [ F"{v+tu)dt, K2{z;u) ^ [ {l-t)F"' {v+tu)dt , Kz{z-u) = [ {l-t)F" {v+tu)dt . 



Thus we can decompose equation (5.1) into linear and nonlinear contributions: 



{dttU - 2sdtzu + s^dzzu) + (l - e^F'{v)) {dtU - sdzu) - dzzU + {se^F"{v)dzV - F' {v)) u 

-2i>'/„.„,\„,/'a„, „S)„,^ I I T^^ I „.\ „^2j^,'„.„,nc),,\„,2 



■Ki{z;u)u{dtu ~ sdzu) + lyK^{z]u) - se^ K2{z;u)dzi') . (5.1) 



Testing equation (5.4 1 against dtu — sdzU yields our first energy estimate (hyperbolic energy): 

\dtu- sdzu\^ dz+\ f \dzufdz+l f {se^F"{iy)dziy- F'{iy))\u\^dz 

+ / {l-e'F'i,.)) \dtu-sdzu\' [ dz{.ss^F"{iy)dz,y-F'{i.))\u\'dz 
Jr ^ JR 

-? / K\{z\u)u\dtu — sdzuf' dz / (K3{z;u) ~ se^K2{z;u)dzi') {dtU ~ sdzu) dz . (5.2) 
Jr Jr 

We are lacking coercivity with respect to norm in the energy dissipation. Testing equation 



dt 1 2 



= £ 



(5.4) against u yields our second energy estimate (parabolic energy): 

' / u{dtu-.sdzu) dz+l f {l-s'^F'{v)) \u\^ dz\ 
Jr 2 Jr J 



d 
dt 



J \dtu - sdzu\^ dz + J \dzu\^dz + J (^^-F"{i^)dzi'- F'iiy)^ \u\^ dz 
= £M Ki{z;u)u'^ {dtu~ sdzu) dz+ / {K3{z;u) - se^K2{z;u)dziy) dz . (5.3) 
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We introduce the following notations for the two energy contributions and the respective quadratic 
dissipations (5.2|,(5.3|: 

I \dtU-sdM^ dz+ \ [ \dM^dz + l [ {se^F"{v)d,v-F'{v)) \u\^ dz , 



E^{t) =e^ I u (dtu - sd,u) dz + - j (l - e'^F'{iy)) \u\'^ dz , 



Qiit) 



(1 - e^F\,y)) \dtu - sdM' + ^ L 9z {se^F"{v)d,v - F'{y)) jup dz , 



QW) = -e' ^ \dtu - sd,\^ + ^ \dM^ + ^ {^-^F"{v)d,v - F' {u)^ \u\^ dz . 

The delicate issue is to control the zeroth-order terms. In particular we define the weights 

A,{z) = se^F"{u)d,v-F'{v), 

9 

A2{z)^^F"{v)d,v~F'{i^). 

They change sign over M. More precisely we have 

lim Ai{z)= lim yl2(z) = > , 
lim Ai{z)= lim ^2(2) = -F'(0) < . 



>+oo 



To circumvent this issue we introduce w{t, z) — e'^^^'^uit, z) as in [2T] and the previous Section 
[Ij where (j){z) is a weight to be determined later (5.9). The new function w{t,z) satisfies the 
following equation: 

e'^dttw - 2e^sdtzw + {2e'^sd,<j> + 1 - e^F'{iy)) dtw + (-.s (l - e^F'{v)) - 2{e^s^ - l)d,<j>) d,w 
+ (e2s2 _ i)d^^^ + (s (1 - e'F^v)) d,cj> + (e^s^ _ ^_q^^^ ^ ^q^^^i^ ^ £''sF"{v)d,v - F' {v)) w 
— e'^Ki{z; u)u {dtw — sdzw) + (K3{z] u) — e'^sK2{z; u)dzi' + e'^sKi{z; u)dz(f>) uw . (5.4) 
We denote the prefactors of dtW, d^w and w as A3, and A5 respectively: 
A3{z) = 2e^sdz<l>+l-e^F\iy), 

Ai{z) = -s (1 - e^F'iiy)) - 2{e^s^ - l)dz(f> , (5.5) 
A,{z) = s (1 - e^F\,^)) dzcp + {e^s^ - 1) + + £\F"{v)dzV - F'{v) . 



Testing (5.4) against dtw, we obtain our third energy estimate: 
d (e^ 



dt 2 



\dtwYdz- 



1 - e^s^ 



\dzw\'^dz+ \ I A^{z)\w\^ dz 



+ / Aj,{z)\dt'w\'^ dz + i Ai{z)dtwdzW dz 
Ki{z;u)u (\dtw\'^ ~ sdtwdzw'j dz+ / [K3{z;u) — sK2{z;u)dzV + sKi{z]u)dz(f)) uwdtw dz , 



Testing (5.4) against w we obtain our last energy estimate: 



d 
di 



wdtuidz+- I A3{z)\w\'^ dz 



\dtw\^dz+{l-e^s'') / \dzw\^ dz + 2se'' / dtwdzW dz + / [ Ar,{z) 



dzAi{z) 



\w\^dz 



[K3{z\u) — sK2{z]u)dzi' + sKi{z\u)dz(t)) uw^ dz+ / Ki[z]u)uw {dfW ~ sdzw) dz , 
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We introduce again useful notations for the two energy contributions and the associated quadratic 
dissipations: 

ETit) = ^ ^ \dtw\' dz + ^ \dM' '^'+IJ^ Mz)\w\^ dz , (5.6) 

E^'{t)^e^ [ wdtwdz+l [ A3{z)\w\^ dz , (5.7) 
Q5"(t)= [ A3{z)\dtw\^ dz + I Ai{z)dtwd,w dz . 



Q2{t)^-£ / \dtw\'dz + {l-e's') / \dM dz + 2e's / dtwd,wdz+ / [ A^iz) 



2 

(5.8) 



To determine (/)(z) we examinate (5.6)-(5.7). We first require the natural condition dz4>{z) > 0. 
This clearly ensures ^3(z) > 1 — e^-F'(O). We examinate the condition ^15(2;) — ^dzA4^{z) > (5.81 
in order to fully determine the weight 4'{z): 

ArAz) - = ,s (1 - s'F'ii.)) + {s\s' - 1) (-a.,</. + + s'sF"{i.)d.i. - F'{v) 

-]^e\sF"{v)d,v+{e^s^~l)d,,<j, 

= {e^s^ - 1) + s (1 - e^F'{v)) d^cj^ + ]^e\sF" {v)d,v - F' {v) . 

This is a second-order equation in the variable dz4>- Maximization of this quantity is achieved 
when 

s{l-e^F'(v)) , , 



We notice that this is equivalent to setting A^{z) = (5.5). Then we obtain 

s'^ (1- e'^F'(v))^ 1 
A,iz) + \e'sF"{.)d.. - F'i.) 

- 4(1 (1 + ^'F'i-f ' 4F'(.)) + \ehF"{.)d.. 

= ^^_^lp,(^o))2 (4^'(0) (1 + ^'F' - 4F'(^) (l + £^^^'(0))') + ^e'sF" i,.)d,i. 
= (^^_Jp,^Qyf (^'(0) - ^'('^)) (1 - ^'^'(O)^'(^^)) + ls'sF"i,.)dz:. . 

We check that As (z) >0 since Vz F'{iy{z)) < F' (0) , F' (0) < l,Vz F"{iy{z)) < and Vz 9^1/(2) < 
0. 

We recall that the exponential decay of ly at +00 is given by the eigenvalue A > 0, where 

s{l-s^F'{0)) 
^ 2(1 -£2^2) • 

Therefore we can rewrite 

Remark 15. As far as we are concerned with linear stability, the energies El" and E'2 contain 
enough information. However proving nonlinear stability requires an additional control of u in L°° 
which can be obtained using i5" and E2 
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2- Combination of the energy estimates. We first examinate the energies and £'2 • We 
clearly have 

EUt) > - i_,2F'(0) 11^*" - "^-"ll' ^^Il"ll2 + ^^MU 

^ II « T ii2 , l-£^f'(0) „ |, 
- ^ l~e^F'{0) " ~ ^ """2 + 1 Il"ll2 

We set 

^ 1 - e^F'iO) 

We have on the one hand 

E^{t) + 5E'2'it)>l\\d,u\\ldz+ [ AG{z)\u\^dz, 

where Aq{z) is defined as 

Mz) = 2^i(^) + ^ — • 

We have on the other hand, 

QUt) + SQlit) > ^ ~ 'T^^^ \\dtu - sd,u\\l + 5\\dM\l + I Aj{z)\u\^ dz , 



where At{z) is defined as 

Mz) = {se'F"ii,)d,,, - F'{u)) + SA^iz) . 

We have both lini^_i._oo ^6(2:) > and limz^^ao At{z) > 0. Accordingly there exists a > 
and zq G M such that 

Vz < zq mm{AQ{z), Aj{z)) > a . 

In order to control the zeroth-order terms over {zo,+oo) we shall use the last two energy 
estimates. First we observe that > zq \u{z)\ = \e-'^'-^^w{z)\ < e--^^^") |w(^)l 

since 4> is increasing. 

We set (j){zo) — without loss of generality. This determines completely together with the 
condition (5.9 1. Wc have 

EUt) + SE^{t)>^\\d,u\\ldz + a f \u\^dz-\\Aee-^H,>,^\\^ f \w\^ dz , 

^ J z<zo J z>zo 

^2 T7l( 



Aje~^^_^ 



. -z>zo 



QUt) + SQ^{t) > ^ ^/'^"^ \\dtu - sd,u\\l + S\\d,u\\l +a f \u\^ dz 

^ J Z<ZQ 

Lemma 16. We have — G L°°(zo,+oo) and — — € L°°{zq, +00). 

A5 Aq 

Proof. The first claim is clearly a consequence of the second claim since Aye""^ G L°°{zq,+oo). 
First we have 

F'{0)-F'{,y) > (^mi^F"^i'. 

Consequently A^ > O^v). Second we recall dziy + Aj^ > (Lemma |9|, so that Vz > zq v{z) > 
O (e""^^). Finally we have Vz > A, thus 

Vz > zo e-^ < p-^(^-^°) <Oiiy)<0 (^5) ■ 



A5{z)\w\^dz. 

00 z>zo 
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We now focus on the second series of energy estimates. We clearly have 



2 1 - £^^'(0) 



1 - £^F'{Q) 



M\2 



> 



2 , l-e^F'{0)^ 



^3^^M2 + 4 Hl2, 



and 



Q^{t) > -e^dMl + (1 - e^s^)\\dM\l - ^^^^\\dtw\\l - ^—^\\dM\l + [ A,{z)\w\^dz 



1 + f2<,2 1-^2,2 



1 - £^S 



?dz. 



oo ^ Z>Zo 



We set 



2£2 



1 + £2s2 



We have on the one hand 
where As{z) is defined as 



1 p2„2 



/- 

Jr 



\\d,w\\2 dz+ / A8{z)\w\'dz 



M^) = ^Mz) + s' > S' 



We have on the other hand, 



A, 



h\w\'^ dz . 



CO '^Z>Zo 



Combining all these estimates, we define E{t) = E^(t) + 5'E^{t) + 5" {E]^{t) + 5E^{t)) and 
Q{t) = Qf{t) + 5'Q^{t) + 5" {Q\{t)+5Q^{t)), where 5" > is defined such as the following 
condition holds true 



5" < y min L^ZM \\A,e-'H^^^,\C 



Are' 



-20 



'-z>zo 



finally obtain our main estimate, 

^E{t)+Q{t)<o( I \u\\dtu-sd,u\''dz+ [ lufdz" 
dt \Jm. Jr j 

+ 0[ [ e-'^\w\\dtw\^ dz + [ e-'''\w\\d,w\'^ dz + [ e-'''\w\^ dz) , (5.11) 
\Js. Jr Jr J 



where 



Eit) > o \\dM\i + / \< dz + \\dM\i + \Mi , 

\ J Z<Zo / 

Q{t) > O (\\dt - sd,u\\l + \\d,u\\l + ( \u\^ dz + \\dM\l + \\dzw\\l + / 

\ J Z<Zq Jz 



^dz] . 



z>zo 
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3- Control of the nonlinear contributions. Our goal is to control the size of the perturbation 
u in L°° . For this purpose we use the embeddings of iJ^(M) into L°°(IR) and 

\\u\\oo<C\\u\\\'^\\dM\l'\ 

\\u\u<c\\u\\l"\\dM\\"- 



We examinate successively the nonlinear contributions. We recall u = e '^w. First we have 



— sdzu\'^ dz < 



\dt - sd-. 



<0[E{tY/^Q{t) 



and similar estimates can be derived for all the contributions in the r.h.s. of (5.11) except for the 
last one. Second we have 



'dz = 



IwP dz - 



Z<Zq 



Z>Zq 



< \H\ 



L°°(-00,2o)II'"IIL2(_oO,2o) 



-0/2, 
|2 



le w\\lHzo,+oo)\Ml^zo,+oo) 



■C\\e- 



'A/2y;l|3/2 



L2(20,-|-00) 



dz 



1/2 

L^iz(}-+oo) 



\w\\lHzo,+oc.) 



We have 



w\ 



L2(_oo,Zo) 



< II"IIl2(-00,2o) ' 



/2 ||2 



z>zo 



< 2 



z>za 



'\dzw\ 



\w\Uz < O 



;\dzcl>\' 



dz <2 



\dM^dz + C 



z>Zo 



Z>Zo 



(E{tY/^Q{t) 



Consequently we obtain 
Finally we get 



dt 

This estimate ensures that the energy is nonincreasing provided that it is initially small enough. 
Indeed there exists a constant C such that jiE{t) + Q{t) < CE{t)'^/'^Q{t). We set c = C-'^/2. If 
initially E'^ < c then the previous differential inequality guarantees that E(t) is decaying and re- 
mains below the level c. Therefore E{t) is positive decaying, and the dissipation Q{t) is integrable. 
This concludes the proof of Theorem [T3| □ 



E{t) + Q{t)<o(E{ty/^Q{t) 



' dz. 
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